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Abstract. We present first results of three-dimensional 
(3D-) calculations of turbulent and degenerate hydrogen- 
burning on top of a C+0 white dwarf of 1.0 M Q . The sim- 
ulations are carried out by means of a code which solves 
Euler's equation for an arbitrary equation of state together 
with a nuclear reaction network and the energy input from 
nuclear reactions on a Cartesian grid covering a fraction of 
the white dwarf's surface and accreted atmosphere. The 
flow patterns we obtain are very different from those of 
earlier 2D simulations using the same initial conditions 
and the same numerical resolution. The possibility of self- 
enrichment of the accreted hydrogen-rich atmosphere with 
carbon and oxygen from the surface layers of the white 
dwarf during the violent phase of the burning is investi- 
gated, and it is demonstrated that self-enrichment pro- 
ceeds too slowly if the accreted gas has near-solar CNO- 
abundances at the onset of the thermonuclear runaway. As 
a result, we do not find a fast nova outburst. This conclu- 
sion remains valid if the initial metallicity of the accreted 
gas is raised by a factor of five. Therefore we conclude that 
fast nova outbursts indeed require huge enrichments of C 
and O, as postulated from spherically symmetric models, 
and that the mechanism which leads to such enhancements 
must operate prior to the outburst. 

Key words: Stars: novae, cataclysmic variables - white 
dwarfs; Physical data and processes: convection - hydro- 
dynamics - nuclear reactions, nucleosynthesis, abundances 



1. Introduction 

The standard model for the outburst of a classical nova 
is a thermonuclear runaway (TNR) in the accreted hydro- 
gen-rich envelope o n top of a white d warf in a close binary 
system ( |5tarrfield| 1989, 1993, 1995; |Truraij 1982, 1990). 
Spherically symmetric models of the runaway, based on 
realistic nuclear reaction rate networks and mixing-length 
theory (MLT) of convection, have been investigated by 



many authors (Starrfield et al. 1974, 1985; Prialnik ct al. 
1978; Mac Donald 1980), and the results were in good 
agreement with observational data, such as the total amount 
of energy released and the metallicity and the abundances 
of the expelled envelope material, provided considerable 
enrichment of the atmosphere with CNO-elements was as- 
sumed. Although the actual numbers depend on the mass 
of the white dwarf and the accretion rate, typical enhance- 
ments of up to a factor of 10 relative to solar values were 
required for models of fast novae. 

Despite of the success of these models, two questions 
remain to be answered: 

— 1. Where does the enrichment come from? 

— 2. What is changed if one avoids the mixing- length 
theory of convection? 

In this paper we concentrate on the second question and 
will eliminate one suggested answer to the first one in 
passing. We shall mainly follow the arguments given in a 
recent work by us (Kercek, Hillebrandt & Truran 1998; 
henceforth referred to as KHT), where these questions 
were addressed by direct two-dimensional simulations. For 
convenience we first summarize here the numerical method 
as well as the main results obtained in KHT. 

In KHT (as well as here) we performed numerical simu- 
lations based upon the hydro-code PROMETHEUS (Fryx- 
ell et al. 1989), a PPM-type code with nuclear reactions 
included. Curvature effects were ignored and the surface 
layers of the white dwarf as well as its envelope were rep- 
resented by a plane-parallel sheet. The advantage of this 
approach was that periodic boundary conditions could 
be used, thereby avoiding common problems of numeri- 
cal simulations of free convection, namely that reflecting 
boundaries may act like a containment and may affect the 
flow patterns in an unphysical way. Effects coming from 
the finite numerical resolution were investigated and it was 
found that, as far as the general properties of the simula- 
tions were concerned, a grid of moderate resolution, e.g. 
220 x 100 (horizontal x vertical) grid points, covering a 
total of 1800km x 1000km of the white dwarf's surface 
and atmosphere, was sufficient to obtain stable results for 
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integral quantities such as the energy generation rate and 
the rate of mixing of carbon and oxygen from the white 
dwarf into the atmosphere. 

In our version of PROMETHEUS, nuclear reactions 
are incorporated by solving together with the hydrody- 
namics a nuclear reaction network including 12 nuclear 
species, i.e., 1 H, 4 He, 12 C, 13 C, 13 N, 14 N, 15 N, 14 0, 15 0, 
16 0, 17 0, and 17 F, linked by reactions described in Wal- 
lace & Woosley (1981). The reaction rates are taken from 
Thielemann (private communication). Following Miiller 
(1986), we solve the network equations and the energy 
source equation simultaneously to avoid numerical insta- 
bilities. 

Since the accretion process prior to the TNR needs 
several 10 5 years in order to get an initial model for the 
multi-dimensional simulations, one has to calculate the ac- 
cretion phase and the slow stages of the burning by means 
of one-dimensional implicit hydro codes. Here as well as in 
KHT we used a model computed to the onset of the vehe- 
ment phase of the TNR by plasner et al. (1997). During 
this early phase, hydrogen-rich matter of solar composi- 
tion (Z = 0.015) is accreted onto the surface of a 1 M Q 
C+O white dwarf at an accretion rate of 5.0-10 _9 M Q yr _1 . 
As the temperature rises at the bottom of the envelope 
convection sets in. Convection is taken into account by the 
MLT in the one dimensional model. The model provides 
the density, temperature, entropy, and composition pro- 
files at the point at which the temperature at the bottom 
of the envelope reaches 10 8 K. The mass of the hydrogen 
shell is then about 2 • 1O _5 M0 and the metallicity is now 
Z = 0.02. This model was mapped onto a radial row of 
the 2D grid in KHT, relaxed for several hundred dynam- 
ical time scales, and then mapped onto the full 2D grid. 
In the present work, the same procedure is applied to a 
3D grid. Finally, an initial perturbation is superimposed 
upon this hydrostatic configuration. 

The main features of the 2D flow fields which showed 
up in KHT were the appearance of small persistent coher- 
ent structures of very high vorticity (and velocity) com- 
pared to the background flow. Although they were not 
fully resolved they caused effects that were essentially in- 
dependent of the resolution. During the early phase of the 
thermonuclear runaway they dominated the flow patterns 
and resulted in very little overshoot and mixing. At late 
times, after steady slow mixing and with increasing nu- 
clear energy production, they became weak, but showed 
up again after hydrogen had mainly been burnt and the 
energy generation rate dropped. The net effect was that 
KHT did find some self-enrichment, but on time-scales 
much longer than in previous calculations. Moreover, for 
initially solar composition of the accreted gas the rise time 
of the temperature was very long, of the order of 1000s, 
and peak temperatures never exceeded 2-10 8 K. Therefore 
these models did not resemble the properties of a fast nova. 

Of course, the question arises if these conclusions would 
remain valid if the assumption of 2D symmetry were drop- 



ped. First, one does not expect to find the small-scale sta- 
ble vortices in 3D because they are a property of plane- 
parallel 2D flows, as was discussed in KHT. On the other 
hand, one also expects the dominant linear scale of convec- 
tive eddies to be smaller in 3D than in 2D and, therefore, 
mixing and convective overshooting could be even less ef- 
ficient than was found by KHT. 

In the following sections we will answer this question. 
In Sect. 2 we briefly describe the parameters of the mod- 
els we have computed. In Sect. 3 we present and discuss 
the results of two sets of computations, one with solar ini- 
tial composition of the accreted gas and a second with 
5 times higher metal abundances. Some conclusions and 
an outlook follow in Sect. 4. We will concentrate on the 
astrophysical aspects of our simulations and leave results 
and questions related to the physics of turbulent combus- 
tion in the so-called well stirred regime to a forthcoming 
paper. 



2. The models 

The computational grid of our present computations cov- 
ers a fraction of 1000km (vertical) x 1800km x 1800km 
(lateral) of the white dwarf's surface and atmosphere by 
100 x 220 x 220 grid points in a Cartesian mesh. This 
is equivalent to the "low resolution" case of KHT. The 
lateral grid is equidistant whereas in vertical direction we 
use a non-equidistant grid with the highest spatial reso- 
lution (5km) in the white dwarf's surface and the bottom 
layers of its atmosphere. As an initial perturbation we in- 
crease the temperature by 1% in 20 zones right on top of 
the white dwarf. In contrast to KHT, we choose a some- 
what larger perturbed area in order to save computer time. 
This, however, is not a problem since the initial conditions 
are a bit artificial anyhow, and we saw in our previous 2D 
simulations that information about the initial perturba- 
tion is lost after about 15s. 

Since we fix the initial conditions according to the 
model of plasner et al. (1997) the only parameter which 
still can be varied within reasonable limits is the initial 
metallicity of the atmosphere. But, because solving the 
reactive Euler-equations in 3D with realistic equations of 
state is "expensive" (on 512 processors of a CRAY T3E 
one time-step needs about 8s CPU-time, and about 2xl0 5 
time-steps are required per simulation), we cannot per- 
form an extensive parameter study. We decided rather to 
choose solar composition for a first simulation (in order to 
compare the results with our earlier 2D computations) and 
to increase the metallicity to 10% in a second run. This 
second simulation is meant to investigate whether or no 
moderately higher CNO abundances lead to significantly 
more violent nuclear burning, thereby changing a "slow" 
nova to a fast one. Although it is certainly of interest to 
increase C and O by another factor of 2 to 3 in order to 
match the conditions which would give a fast nova in ID 
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simulations, we did not do it because of obvious inconsis- 
tencies of our initial model in that case. 

Combustion in novae is peculiar because convective 
overturn times are much shorter than the nuclear reac- 
tion time-scale. Therefore, in general, stirring proceeds 
much faster then nuclear burning. In our approach to this 
problem, namely by direct numerical simulations, we re- 
solve the large scale inhomogeneities and velocity fluctu- 
ations well, but the small scale dissipation and mixing is 
dominated by numerical discretization effects. This is not 
a problem as far as the energy transport by convective 
motions is concerned because the energy flux is predomi- 
nantly carried by the largest eddies. Mixing by turbulent 
motions, however, happens on very small scales which are 
definitely not resolved in our computations. Therefore we 
will overestimate the mixing due to the numerical diffu- 
sion of matter on the grid scale. But since, as we shall 
show, even with our treatment mixing is too slow for a 
fast nova to occur, direct simulations are a conservative 
approach to the problem under consideration here. 

3. Results of the 3D simulations 

3.1. The case of near-solar metallicity (Z=0.02) 

In this subsection we present and discuss the results we 
have obtained for an atmosphere of nearly solar metallic- 
ity, the main emphasis being on a comparison with KHT. 
The initial abundances were obtained by evolving the ac- 
creted atmosphere of a given metallicity through steady 
hydrogen burning to the vehement stage of the TNR. We 
shall first analyse the evolution of the convective flow 
fields, and then discuss integral quantities such as the lat- 
erally averaged temperature and the nuclear composition. 
We computed this model over about 400s. 

3.1.1. The evolution of the flow fields 

In order to visualize our results we present them in the 
form of 2-dimensional cuts and iso-surfaccs in three di- 
mensions. The 2D-cuts show, colour-coded, the absolute 
values of the velocity as a function of time. They indi- 
cate where in the atmosphere of the white dwarf most of 
the convective transport is happening and can easily be 
compared with our earlier 2D results. In contrast, the iso- 
surfaces of the absolute values of the velocity field at fixed 
times are used to discuss the characteristic structure and 
scales of the convective eddies. 

Figure [I] gives the temporal evolution of the velocities 
for a typical 2D cut. As in the 2D simulations of KHT, 
the surface layer is ignited by an "ignition string" per- 
pendicular to the plane shown. After a couple of seconds 
the entire surface layer is burning and the first convective 
eddies form. This stage is reached a little earlier than in 
the 2D simulations, mainly because a larger volume was 



perturbed initially. From then on the evolution is very 
different from the 2D case. No regular structures appear. 
In contrast, the flow field is very irregular. Violent erup- 
tions occur occasionally at later times, followed by periods 
of rather quiet burning. Towards the end of the compu- 
tations, convective motions extend all the way into the 
upper atmosphere. 

The evolution of the convective patterns can be seen 
more clearly in the iso-surface plots Figs. || to [I| Figure 
||, taken at 50s, shows that most of the fast eddies in the 
beginning are confined to a narrow layer near to the white 
dwarf's surface, extending over roughly 150km. This burn- 
ing layer is well separated from an upper convective zone 
with large but slow eddies. But even the convective mo- 
tions near the surface are slow in comparison to the sound 
velocity which is above 10 8 cm/s there. During this phase, 
the burning matter does not expand significantly and also 
the upper atmosphere is only moderately heated from be- 
low. Consequently, the temperature rises in the burning 
region above 10 8 K, and the burning becomes more vio- 
lent. As is shown in Fig. ||, after about 100s, convective 
velocities begin to exceed 10 7 cm/s and some eddies start 
to penetrate into the upper atmosphere. From then on 
very efficient heat transport by convection leads to an ex- 
pansion of the atmosphere, and energy going into lift-off, 
and adiabatic energy losses balance essentially gains from 
nuclear burning. After 200s, finally, the entire atmosphere 
is well stirred by convective eddies of all sizes with typical 
velocities around 6T0 6 cm/s (see Fig. |^). It is obvious that 
the flow fields in 3D are very different from the ones in 
2D, as computed in KHT, and one expects that this will 
also effect integral quantities. 

When we stopped this set of computations after 400s, 
the nuclear energy generation rate is approximately con- 
stant throughout the white dwarf's atmosphere and a stage 
of nearly steady hydrogen burning is reached. The aver- 
age nuclear energy generation rate has come down to a few 
times 10 13 erg/g/s, missing the conditions for a fast nova 
by far. The total energy released by nuclear reactions is 
5.8T0 44 erg which should be compared with the value of 
2.M0 45 erg for our 2D simulation at the same time. 

3.1.2. Temperature evolution, energy generation and 
convective mixing 

The temperature evolution of this near-solar metallicity 
model is depicted in Figs. |^ and [| and the corresponding 
energy generation rates are given in Fig. 0. 

In Fig. |5| we picked the temperature of the hottest zone 
and plotted it as a function of time. The first short drop 
reflects the ignition phase. The following rapid rise to the 
peak value of about 1.1-10 8 K is due to degenerate nuclear 
burning without efficient convective energy transport and 
expansion, as was discussed in the previous subsection. 
This phase is followed by some expansion and adiabatic 
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cooling, accompanied by the onset of more efficient con- 
vective energy transport and mixing. Consequently, the 
maximum temperature rises again, but now its peak value 
is found a bit further out in the atmosphere. These tem- 
peratures are typically 10 to 20 % below what we have ob- 
tained in our 2D simulations mainly because, as we shall 
see, mixing of fresh C and O from the white dwarf into 
the hydrogen-rich atmosphere is slower in 3D than in 2D. 

Figure ^ shows snapshots of the temperature distribu- 
tion, averaged over horizontal planes. With the exception 
of a clear spike at 98s, indicating nuclear burning in a well 
confined surface layer, the distributions are very smooth, 
reflecting the nearly adiabatic average temperature gradi- 
ent of a fully convective stratified atmosphere. 

These effects are more clearly seen in the vertical pro- 
files of the nuclear energy generation rate, again averaged 
over horizontal planes, given in Fig. |^. For the first about 
100s we observe a rapid increase of the energy genera- 
tion rate in the layers just on top of the white dwarf, and 
convective energy transport is insufficient to induce fast 
burning in the outer parts. Towards the end of this simu- 
lations one can see that a stage of nearly steady burning- 
is reached throughout most of the atmosphere. 

Finally, in Figs. || and ^| we show the mass fraction of 
12 C and the total metallicity for the same snapshots as in 
the previous figures. The sharp discontinuity at 5.5-18 8 cm 
in the initial values presents the white dwarf's surface. 
The initial mass fraction of 12 C is about 0.06 % only 
because steady hydrogen burning prior to the TNR has 
transfered most of the originally solar CNO isotopes into 
14 N, whereas the metallicity is untouched. One can see in 
Fig. H as well as in Fig. || that fresh 12 C is mixed into the 
atmosphere by convective undershoot at a very moderate 
rate only. By the time the electron degeneracy is largely 
removed (after about 100s) mixing has just doubled the 
metallicity in most of the burning zones. This is the main 
reason why we did not obtain a violent outburst. 

3.1.3. Mass loss and nucleosynthesis 

In concluding this section, we want to discuss briefly pos- 
sible observable consequences of this model. Peak out- 
flow velocities are reached after 170s and they amount 
to 5-10 6 cm/s at a radius of 6400km which is far below the 
escape velocity from the white dwarf of about 6-10 8 cm/s 
(see Fig. |o|). Typically, the velocities are also a factor of 
two below those we obtained in the 2D run at the same 
instants of time. Therefore, although we are loosing mass 
from our computational grid by outflow (see Fig. [ll]), there 
is no direct mass ejection from the white dwarf. In con- 
trast, mass loss will happen in form of a wind from the 
outer atmosphere which was not included in our compu- 
tations. 

Conclusions concerning the chemical composition of 
the wind, therefore, are rather uncertain. Nevertheless we 



can obtain rough estimates by looking at the abundances 
we find in the computational domain. As can be seen in 
Figs, [l^ through [l^, after about 400s (in reality already af- 
ter 200s) the composition of freshly synthesized (radioac- 
tive) isotopes is homogeneously distributed throughout 
the white dwarf's atmosphere. Since we expect that fast 
convective mixing will proceed all the way out to optically 
thin regions one can also expect that this composition is 
a fair representation of the matter which is lost from the 
white dwarf ejected in the wind. At late times, radioactive 
decay heating by these isotopes dominates the energy pro- 
duction in all but the innermost layers of the atmosphere, 
but the overall enrichment of CNO isotopes is only moder- 
ate, probably not more than a factor of five to six relative 
to the initial values. 



3.2. A model with Z=0.1 

In one-dimensional nova simulations initial metal abun- 
dances largely in excess of their solar values are required 
to drive a fast outburst. Since we have demonstrated that 
self-enrichment of the atmosphere by convective under- 
shoot and dredge-up of matter from the white dwarf does 
not happen during the outburst, we decided to run a 
model in which the atmosphere was already enriched from 
the start. In order to avoid large inconsistencies (this model 
was not evolved in the same way as the low-metallicity 
one), we increased Z to 0.1 only and constructed a new 
equilibrium model which had the same temperature struc- 
ture as before. Burning was then ignited by raising the 
temperature by 1% in 20 zones of the bottom layer of 
the accreted atmosphere, and all other properties of the 
model were identical to the previous one. Of course, an 
enrichment by a factor of five is still below what is com- 
monly assumed in ID simulations but one might hope for 
more violent burning and, therefore, also faster mixing of 
C and O already in this case, if convective motions are 
treated properly. However, as we shall demonstrate, this 
is not the case. We find that the evolution proceeds con- 
siderably faster than in the low-metallicity case, but the 
final outcome is not too different in both cases. 

Since the flow patterns look very similar to those shown 
in Fig. 0, with the exception that already after about 80s 
velocities above 10 7 cm/s are observed, and after about 
150s peak velocities reach 2-10 7 cm/s, nearly a factor of 2 
higher than in the low-metallicity run. The second effect is 
that peak temperatures are now a bit higher (see Figs. [i~5| 
and [if], to be compared with Figs, [sj and ||), and also the 
energy generation rate exceeds the previous one, but not 
by much. After 160s a state of burning is reached which is 
nearly indistinguishable from the low-metallicity case at 
400 seconds. Even the metallicity in the atmosphere (Fig. 
|l7| ) is only moderately higher than before. 

We therefore conclude that, although originally higher 
CNO-abundances speed up the thermonuclear burning, 
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going to 5 times solar metallicity is still insufficient for 
a fast nova. The reason seems to be that, despite the fact 
that the consumption of hydrogen is faster in the begin- 
ning, the time scale for the temperature rise is still too 
long, i.e. on the order of 100s. This, in turn, means that, 
as before, the atmosphere begins to expand and cool once 
the temperature exceeds 10 8 K, much too low for a fast 
nova. Very fast mixing of white dwarf material into the 
atmosphere during the early phase of the TNR could, in 
principle, overcome this problem but we do not find such 
a mixing. In contrast, the 3D simulations mix even less 
than the 2D models of KHT. 



4. Summary and conclusions 

Wc have presented the first 3-dimensional simulations of 
thermonuclear models of classical novae. We have shown 
that the results differ, in several aspects, considerably 
from those of previous 2-dimcnsional simulations and we 
conclude that this and similar problems have to be carried 
out in 3D in order to give reliable results. The main reason 
is that the flow fields in 3D are very different from those 
in 2D and, therefore, all quantities which depend on them, 
such as convective energy transport and mixing, will also 
differ in both cases. For example, in 3D we find, not un- 
expectedly, more power in small scale motions than in our 
previous 2D simulations which leads to less convective un- 
dershoot into the white dwarf's surface and consequently 
less dredge-up of C and O. This, in turn, makes explosive 
hydrogen burning less violent, weakening the chances for 
getting a fast nova that way. 

Therefore we arrive at the conclusion that self-enrich- 
ment of the accreted atmosphere with C and O during the 
outburst is very unlikely, if not impossible, and our simu- 
lations rule out one of the suggestions for changing a slow 
nova into a fast one. A more likely solution to this prob- 
lem seems to be large enrichment prior to the outburst by 
either shear-induced instabilities at the interface between 
the white dwarf and its atmosphere (like in an accretion 
belt) or some kind of diffusive or convective mixing during 
the long quiet accretion phase. 

Finally, we want to point out that the type of numer- 
ical simulations we have performed leave room for im- 
provements. Firstly, as was mentioned already in the intro- 
duction, in our numerical method mixing on small scales 
is due to numerical diffusion. It is difficult to estimate 
quantitatively its effect. It appears to be safer to include 
mixing by small scale turbulence explicitly by a subgrid 
model, as was done for the propagation of nuclear flames 
by Niemeyer & Hillcbrandt (1995). Secondly, the kind of 
problems wc have tackled here with brute force, namely a 
problem in which the typical velocities are far below the 
sound velocity, should be better approached by means of 
(nearly incompressible) implicit hydrodynamics schemes 
which are presently being developed. It is obvious that 



one would like to explore a larger fraction of the param- 
eter space, but this can only be done with more efficient 
codes. 
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Fig. 1. Snapshots of the velocity field of the three-dimensional run with an initial metallicity of Z = 0.02. Shown are 
vertical slices through the calculation domain. The absolute values of the velocity at each point are color coded. Note 
that the "radial" coordinate is logarithmic above 6500 km. 



Fig. 2. Iso-surfaces of the absolute value of the velocity for the 3D simulation with an initial metallicity of Z = 0.02 
at t — 50 seconds. The (800km) 3 cube shown in this figure is only a fraction of the total computational domain. 
Two distinctly different convective layers can be seen: one of high velocity near to the surface of the white dwarf, and 
a second one of low velocities in the outer atmosphere. 



Fig. 3. Same as Fig. 2, but at t — 100 seconds. It can be seen that typical convective velocities still decrease with 
increasing distance from the white dwarf's surface, but the different layers begin to mix. Those close to the surface 
are heavily stirred by nuclear energy generation. 



Fig. 4. Same as Fig. 2, but at t — 200 seconds. The convective eddies no longer show a clear structure. The entire 
computational domain is well stirred and covered by eddies of all sizes and different velocities. 
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Fig. 5. Time evolution of the horizontally averaged temperature in the hottest layer of the atmosphere for the 3D 
run with Z = 0.02. 




5.6x10 u 5.8x10" 6.0x10 u 

Radius (cm) 



Fig. 6. Horizontally averaged vertical temperature profiles for the 3D run with Z = 0.02 at s (solid line), 98 s 
(dashed-dotted line), and 400 s (dashed line), respectively. 
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Fig. 7. Same as Fig. ||, but for the horizontally averaged vertical profile of the energy generation rate. 
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Fig. 8. Same as Fig. 6, but for the horizontally averaged vertical profiles of the 12 C mass fraction. 
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Fig. 9. Same as Fig. 6, but for the horizontally averaged vertical profiles of the mctallicity. 
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Fig. 10. Horizontally averaged vertical profiles of the vertical component of the velocity at four different times for 
the 3D simulation with Z = 0.02. 
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Fig. 11. Horizontally averaged vertical profiles of the density for the run with Z = 0.02 (0 s solid line, 200 s dashed- 
dottcd line, 400 s dashed line, respectively). 




Fig. 12. Snapshots of the horizontally averaged vertical profiles of the 15 mass fraction for the case Z = 0.02. 
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Fig. 13. Same as Fig. ||, but for the 14 mass fraction. 
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Fig. 14. Same as Fig. ||, but for the 17 F mass fraction. 
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Fig. 15. Temperature in the hottest shell of the atmosphere for the 3D run with an initial metallicity of Z = 0.1. 
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Fig. 16. Snapshots of the horizontally averaged vertical profiles of the temperature for the run with Z = 0.1 initially, 
at s (solid line), 78 s (dashed-dotted line) and 160 s (dashed line). 
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Fig. 17. Same as Fig. |6[ but for the metallicity. 
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